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SECTION 1 
INTRODUCTION 


Ground wave propagation is possible at all radio frequencies but 
often neglected in the HF band (3-30 MHz) because of the much longer paths 
possible with sky wave propagation. Typical ground wave propagation paths 
are measured in tens to hundreds of kilometers while typical sky wave 
paths extend to thousands of kilometers. The purpose of this paper is to 
consider various aspects of HF ground wave propagation over smooth and 
irregular terrain. Although HF frequencies are emphasized, many of the 
results are applicable to a much wider range of frequencies. 

A military situation involving a large nuclear exchange is one 
example where ground waves may be considerably more effective than sky 
waves. Nuclear explosions cause severe D-region absorption which can 
reduce or eliminate conventional HF sky wave modes. These explosions are 
expected to have little or no effect upon ground waves. While sky waves 
can be degraded by attacking the propagation medium (the ionosphere) such 
a tactic is nearly impossible with ground waves. Ground waves propagate 
along the ground/air interface and are only sensitive to substantial 
changes in this interface. Even a large nuclear exchange would not be 
expected to greatly alter the basic topography of a region and hence, the 
ground wave propagation. 

On the other hand, ground wave propagation is vulnerable to jam¬ 
ming or to direct destruction of the terminals as is sky wave propagation. 
However, jamming of ground wave circuits can be made difficult by taking 
advantage of the limited range of ground waves. If sky wave propagation 
is suppressed by proper choice of frequencies and antennas, then an enemy 
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must place a jammer near or within sight of the terminals being used to 
effectively jam them. 

Ground wave propagation generally requires as much or more trans¬ 
mitter power and as good or better antennas than sky wave propagation. 
Higher transmitter power enhances enemy direction finding and eavesdrop¬ 
ping but is compensated by the relatively small ground wave coverage 
regions. Because ground wave propagation is better at lower frequencies, 
it is tempting to consider operation at the low end of the HF band or even 
into the MF. However, efficient antennas at these frequencies are large 
and probably impractical for certain mobile applications. For these 
applications, it may be better to operate at middle HF frequencies where 
efficient antennas of moderate size are possible. 

In general, ground wave propagation offers worthwhile capabili¬ 
ties for short distance military communications systems intended to sur¬ 
vive a nuclear attack. Network communications systems with many closely 
spaced nodes may find ground waves attractive. 

Because the subject of ground wave propagation is broad and 
complex, this report concentrates on those aspects which have recently 
been studied for DNA. An extensive survey of the ground wave literature 
has been made leading to the selection of a theoretical formalism appro¬ 
priate to computer coding. The formalism chosen was developed by Hendri- 
fcus Bremmer over 30 years ago and applies to propagation over a smooth 
spherical earth. This formalism has been coded in Fortran IV and is 
reported here. The computer routine has been kept relatively simple and 
does not yet include corrections for atmospheric refraction or for rough 
earth cases. These subjects are, however, discussed in some detail. 
Finally, the lack of experimental data is discussed along with a sugges¬ 
tion for a future experiment. 

This work was supported by two DNA contracts: DNA 001-80-C-0022 
and DNA 001-80-C-0225. 
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SECTION 2 

GENERAL CONSIDERATIONS 


There is often confusion as to the physical reality of ground 
waves (or surface or terrestrial waves). They are that portion of elec¬ 
tromagnetic waves propagated through space and affected by the presence of 
the ground. They do not include any portion of the wave reflected from 
anything other than the ground. For example, ionospheric sky waves or 
tropospheric waves are separate phenomena. Line-of-sight propagation 
between two spacecraft in empty space away from the earth is also a separ¬ 
ate phenomenon although "line-of-sight" propagation between two aircraft 
aloft is technically ground wave propagation. As the aircraft fly higher, 
ground wave propagation approaches line-of-sight propagation. In other 
words, the dielectric properties (dielectric constant and conductivity) as 
well as the topology of the ground influence the propagation of electro¬ 
magnetic waves "near" the ground. "Near" is defined in terms of the 
number of wavelengths that the transmitter and receiver are above the 
ground. At two meters above the ground, visible light is hardly affected 
at all by the ground while HF transmissions are strongly affected. 

Ground wave propagation is a diffraction effect but not one of 
the simpler diffraction effects often considered in the study of optics.2 
Although the general concept of diffraction around an edge applies to 
ground wave propagation around the earth, the mathematical formalism is 
considerably more complex for ground wave propagation. At one time, this 
subject was considered to be among the most difficult in theoretical 
physics. To obtain any solution at all requires several assumptions. 
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Arnold Sommerfeld solved the general problem of radio propaga¬ 
tion from a short vertical antenna over a finitely conducting plane earth 
in 1909.3 The rigorous results obtained by Sommerfeld were applicable 
to short transmission distances where the curvature of the earth could be 
ignored. It was already known that the solution for a spherical earth 
could be obtained in terms of a series of zonal harmonics. Unfortunately, 
this series converged so poorly that no practical results could be obtain¬ 
ed for radio waves. 

In 1918, Watson^ developed a transformation with the aid of an 
integral in the complex plane which transformed the solution in terms of 
zonal harmonics into a residue series which converged rapidly and could be 
evaluated numerically. From 1926 through 1931, Sorrcnerfeld, Van der Pol, 
Niessen, and Wise obtained several independent solutions of the problem 
using Watson's transformation. During this period, a fundamental error 
was also uncovered and corrected in Sommerfeld's original work. In the 
mid-1930's, Bremmer and Van der Pol in Europe and Norton in the United 
States extended the early work to the point where it became practical to 
use. 5 »6,7,8,9,10 They simplified the mathematics, performed extensive 
numerical computations, and helped to better define the physical nature of 
ground waves. After the Second World War, many other investigators 
considered various aspects of the problem such as an inhomogeneous or 
rough earth. The most generally useful results, however, are still those 
of Bremmer or Norton. We have arbitrarily chosen to follow Bremmer's 
formalism^ because it appears to be the easiest to use. 

Bremmer's calculations reveal several general characteristics of 
ground waves. For instance, they propagate much better over water than 
over land, and better over wet land than dry land because the high conduc¬ 
tivity of water and wet land produces less absorption. Ground waves are 
more influenced by the conductivity of the ground (soil type) than by the 
curvature of the earth (for transmitter and receiver on the earth). They 
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are not abruptly changed by the horizon. As with most diffraction 
processes, there are no distinct shadows behind obstacles such as the 
horizon. Also, long wavelengths propagate farther around the earth than 
do short wavelengths. Ground waves are usually vertically polarized 
because horizontal polarization is readily short-circuited by the earth. 
Hence, horizontally polarized antennas such as might be used for sky wave 
propagation are relatively ineffective for ground wave propagation. 








SECTION 3 

BREMMER'S EQUATIONS/GROUND WAVE PROGRAM 


This section presents a seminary of Bremmer's approach to ground 
wave theory and shows how this theory has been used to create a simple 
computer routine for computing ground wave field intensities around a 
smooth spherical earth. 

As mentioned in the last section, the rigorous theory for wave 
propagation around a sphere is a special application of diffraction 
theory. It is mathematically difficult because of the finite size of the 
sphere. Much simpler solutions are possible for a very small or large 
sphere. The general solution in terms of a series of zonal harmonics is 
given by^ 

H ■ ^ [— ♦ ik o hw\ 7TT7rV )(k o b) 4 l) (V) p n< c ° s8 n 

cb s n =0 4 U k o a ) (1) 

where H is the magnitude of Hertzian vector. This Hertzian vector is 
related to the conventional vector potential. A, by the relation 


Hr = _ 
v 


A = _ 


aA 

IT 


( 2 ) 


From a knowledge of the Hertzian vector, the electric and magnetic field 
components follow directly. 
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The effective reflection coefficient in (1) is defined by 


and the functions c n ( L )(x) and f n (x) are given by 


where H and J are the conventional Hankel and Bessel functions. Other 
quantities in (1) are defined as 

IL = electric dipole moment of the assumed small 
vertical dipole antenna 

s = distance between transmitter and receiver 

b = radial distance of antenna from center of assumed 
spherical coordinate system 

k 0 ,k! = radial wave numbers outside and inside the earth 
a = radius of sphere (in this case, the earth). 


The numerical evaluation of (1) is very difficult and for a long 
time remained unknown. The problem is the convergence of the series for a 
finite radius, a. If a is small compared to the wavelength, the series 
converges rapidly, and the first term yields Rayleigh scattering. The 
greater the value of the parameter k x a = 2*a/x, the more slowly the series 
converges. The most significant terms are those of order n - kja. In the 
radio case, k L a varies from 10 3 to 10 9 requiring the summation of a great 
number of terms. 











Because such a summation is impractical, Watson developed a 
transformation which converted (1) into a contour integral in the complex 
n-plane. This reduced the problem to a summation of residues which can be 
mastered numerically. Bremmer introduced approximations for the resulting 
Hankel functions and zonal harmonics leading to a simpler expression for 
the r.m.s. value of the electric field at a distance, D km , from a short 
vertical dipole antenna.1 That expression is shown below with two 
significant modifications to bring these computations in line with the 
conventions of the CCIR (see Appendix). The electric field has been 
normalized to give a value of 346.4 mv/m at 1 km if the assumed dipole 
transmitting antenna is placed on a perfectly conducting plane earth. The 
same dipole antenna will give a field of 173.2 mv/m at 1 km in free space 
which is the free space field from a 1 kw isotropic antenna. Bremmer and 
others use 300 mv/m instead of 346.4 which corresponds to the field from a 
dipole antenna radiating 1 kw over a perfectly conducting plane earth. 
The antenna assumed here radiates 4/3 kw under these conditions. Another 
modification to Bremmer 1 s equations involves the distance D km . While 
Bremmer uses straight line and great circle distance equivalently, the 
following expressions use the latter only. 


346.4-y2irx /\ (hj+a) 2 + (h 2 +a) 2 - 2(h 1 +a) (h 2 +a) cos (D^/a)] 


x I I MM f s (M e lTsX /(2T,-l/6 2 )| mv/m (6) 

s=0 5,1 

where 

D. = distance along the surface of the earth from transmit¬ 
ter to receiver in km. 


x • -053693 0^/^1/3 (7) 

x m = wavelength in m 
f s (hi) = height gain factor 
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x^-2, s \l/2 

, / 


^1/3' [i (*, 2 -2t s ) s / 2 ) 
H,/3 l,) [j C-2t s ] 3/2 ] 


(8) 


Xl 2 = .03674 h,/x2/3 

h : = height of transmitter above ground in m 

h 2 = height of receiver in m 

a - K i(135*-(> e ) 
e e 


K 


e 


.002924 X m 
m 


1/3 * 3.6«10 25 aJ ^ 

\Ti e-l) 2 + 3.6 x10 25 o 2 xj 


i|> = tan -1 (--- 

6x10 12 0 A 

e m 


) - _ tan- 1 (- 

2 6x10 12 


-1 


\n 


The values of x s follow from 

t = Re r + i Im t 
s s s 

and assuming K e is small: 

Re t 0 = .928 + K cos(45* + Y ) + 1.237 K 3 cos(75* + 3 y) 

U G 6 G G 

- .5 K " cos(4Y ) - 2.755 K 5 cos(75° - 5 yJ ... 
e e e e 

Re Tl = 1.622 + K e cos(45° + Y g ) + 2.163 K g 3 cos(75* + 3Y g ) 

- .5 K e 4 cos{4Y e ) - 8.422 K g 5 cos(75* - 5Y e ) ... 


(9) 

( 10 ) 


(ID 


( 12 ) 
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= 2.694 + K e cos (45° + «F e ) + 3.592 K g 3 cos (75° + 3f g ) 

- .5 K g 4 cos (4» e ) - 23.227 K g 5 cos (75° - 5^) ... 

= 1.116 (s + 3/4) 3 + K cos (45° + ? ) 

G 6 

+ 1.488 (s + 3 / 4 ) 2 /3 k.^3 cos (75° + 3¥g) 

- .5 K g 4 cos (4^) 

- 3.987 (s + 3/4)*♦/3 K g 5 cos (75° - 5^) ... for s > 3 

= 1.607 - K g sin (45° + ¥ g ) - 1.237 K g 3 sin (75° + 3* 0 ) 

+ .5 K g 4 cos (4¥ g ) - 2.755 K g 5 cos (75° - 5¥ ) ... 

= 2.810 - K e sin (45° + ¥ g ) - 2.163 K g 3 sin (75° + 3¥ e ) 

+ .5 K g 4 sin (4» g ) - 8.422 K g 5 sin (75° - 5^) ... 

= 3.795 - K g sin (45° + ¥ g ) - 2.921 K g 3 cos (75° + 3^) 

+ .5 K 4 sin (4v ) - 15.36 K 5 sin (75° - 5¥ ) ... 

c e 6 G 

= 4.663 + K g sin (45° + ¥ fi ) + 3.592 K g 3 cos (75° + 3fg) 

+ .5 K e 4 sin (4¥ e ) - 23.227 K g 5 sin (75° - 5¥ g ) ... 
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1 -932(s + 3/4) 2//f3 - K sin (45° + y ) 
e e 

- 1.488(s + 3/4) 2 /3 k 3 s -j n ( 75 ° + 3 ^, \ 

e e 

+ .5 K 4 sin (4'F ) 
e e 

- 3.987(s + 3/4)‘t/ 3 K g 5 sin (75° - 5^) ... for s > 3 


For K e large: 


.4043 + .6183 


cos (15° - r e ) sin(2v p ) 

-!l_ - .2364 -- e . 


cos (15° + 30 cos (60° - 40 

- .0533 -51 + .00226 _._ e 

K e 3 K e " 


1.288 * .194 C ° S (15 ‘ - . , 0073 S" V'*) 


cos (15° - 30 cos (60° - 40 

+ .0120 -11 - .00160 ---il 


Re T s - 


1.116 (s + 1 / 4 )2/3 + _ l 2241 C ° S (15 I 4 ' e ) 

(S + 1/4) 2/3 Ke 


for s > 1 


.7003 - .6183 11 1 15 ' - ^ ♦ .2364 ffilZi 

Kq k 2 


- .0533 


sin (15° + 30 


sin (60° - 40 


1 

1 


15 





2.232 - .1940 


+ .0073 


Im tj 


Im T s 


sin (15° - v) 
6 


sin (15* + 3^) 

+ .0120 __ + . 


K e 3 


00160 


1.932 (s + 1/4) 


2/ 3 


.2241 


(s + 1 / 4 ) 2 /3 


cos (2^) 

" Ke 2 
sin (60* - 4f e ) 

*7 

sin (15* - f e ) 

K e for s>l 


These expressions apply to a vertical dipole transmitting 
antenna and to a separation between transmitter and receiver of at least 


D|<m ^ ^ Am 1 ^ 3 * 

For a horizontal dipole antenna, the above equations can be used 
(in the direction of the maximum field) providing that 6 e , K e , and i|> e are 
replaced by s m , and 90’ - f m . These latter quantities are defined as 
follows: 


Km e 


1(45* + fj 


(13) 


K m = .002924 


y / 3 _ 

*W(e - l) 2 + 3.6 x 10 25 o e 2 Xm 2 


(14) 


1 tan- 1 (_ til _) 

2 6 x 10 lz ^ 


(15) 


It should also be noted that a e is defined in terms of e.m.u. 
following Bremmer's preference. In terms of more customary usage, this 
can be written as 

0 = 10 11 o o 

e n 

e.m.u. mhos/m 
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The height gain factors in (6) and (8) are made up of Hankel 
functions. Using Bremmer's approximation for these functions, they can be 
written as 


a) for h lm > 50 A,,, 2 / 3 (approximately): 

f s (M - je 1 ' ,V4 * J 'S' ' 2 ^ )3/2 ' „ “ 6) 

5 e V - 2t s < 

", _ j .2038 .3342 

U*- 2x s )^' (x t 2 - 2 t $ ) 3 ‘ 


where 


,[W4 - j(x L 2 - ^ T s ) 3/2 ] 


1 + i 


.2083 


( Xi 2 - 2 t s )3/2 


A 0 = .3582 e 1 120' 

A l = .3129 e-i 60' 

A 2 = .2903 e 1 120’ 

A 3 = .2760 e- 1 ' 60# 

A = .3440 -L ±1 _ 

(s + 3/4) 


-in/3 


for s > 3 


(17) 


b) for h < 50 A m 2 ^ 3 (approximately): 


f (h,) = 1 ♦ 6.283 (—J— - I) - 39.48 * e T *> (h)*... 

x‘ /3 S e X X X-/S 5 e X 

where 

= 4xlQ 7 
*m 


(18) 
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The height gain factor for the receiver f s (h 2 ) can be computed 
as above by substituting h 2 for h r . Similarly, the height gain factors 
for a horizontal dipole can be computed from the equations above by 
substituting ^ for 6 e . Bremmer maintains that methods (a) and (b) 
above are sometimes equally suitable. 

Using equations (6) through (18), a ground wave program has been 
constructed and is given in Table 1. By modifying the input and output, 
this program has also been converted into a subroutine for the nuclear 
effects code known as HFNET. The necessary inputs (including units) are 
shown in the program listing. This program applies to ground wave propa¬ 
gation around a smooth homogeneous spherical earth. It does not contain 
any provisions for a rough or inhomogeneous earth. It also does not 
contain the effects of atmospheric refraction. Despite these limitations, 
the program accurately computes ground wave field strengths in many prac¬ 
tical cases where there are no large departures from the assumptions built 
into the program. The program is not limited to the HF band but applies 
to a wide range of radio frequencies. It has been checked out from 15 KHz 
to 100 GHz but may be applicable to an even wider range of frequencies. 
It is also applicable to a wide range of separations between transmitter 
and receiver. As will be shown below, the routine does fail under some 
circumstances where the transmitter and receiver are close together (in 
terms of number of wavelengths) or high above the earth in full view of 
each other. Under these conditions, a "geometric optical" model is more 
appropriate and will be added to the program in the future. Although the 
residue series model used in this program is relatively less efficient 
when the receiver is in view of the transmitter, the results are still 
accurate except as noted. Many more terms in the residue series need to 
be summed »rf>en the receiver is above the transmitter's geometrical hori¬ 
zon. 


Extensive checkout of the ground wave program has been accom¬ 
plished and is shown in Figures 1 through 8. (All of this checkout was 
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Table 1. Ground wave program listing. 

C PROGRAMi GROUND WAVE 


C 

C 

C 

C 

C 

C 

C 

0001 
0002 
0003 
0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

C 

C 

r 

c 

c 

c 


**************************************** 

BREMMER GROUND WAVE MODEL 
CODED BY GORDON J. FOLKS 4/80 

**************************************** 

DIMENSION TAUR<0J3)r TAUI(0J3>r HT(2)r CSQ(2> 

REAL KrK2rK3rK4'K5rLAMBDA 

COMPLEX DELT A r Jr A(013)r DIFr Fir F2r AX 
COMPLEX SUM r F3r SUMXr SSDIFr HGF< 2)r TAU 
OPEN <UNIT“1rNAME='GND.DAT'rTYP£='OLD'rREADONLY > 

OPEN (UNI!=2r NAME”'GND•OUT'r TYPE®'NEW') 

Cl * 1./3. 

C2 - 2. # Cl 
C3 « 4. * Cl 
C4 - Cl / 2. 

J = CMPLX< 0. r 1.) 

PI * 3.1415976535 
RAD a 6.371E3 

**************************************** 

INPUT SECTION 

**************************************** 

READ TRANSMITTER POWER IN KILOWATTSr TRANSMITTER WAVELENGTH 
IN METERS» AND GREAT CIRCLE DISTANCE ALONG THE SURFACE OF 
THE EARTH FROM TRANSMITTER TO RECEIVER IN MlOMETERS 


0014 


READ <lr*> PWRr LAMBDA» D 


0015 


C 

C 

C 


READ HEIGHTS OF TRANSMITTER AND RECEIVER ABOVE THE GROUND 
IN METERS AND WHETHER THEY ARE VERTICALLY OR HORIZONTALLY 
P0LARI7ED <0 - VERTICALr 1 - HORIZONTAL) 

READ <1r *) HT 9 NPOL 


C 

C 

C 

C 

0016 

c 


READ GROUND CHARACTERISTICS 

EPSLN RELATIVE DIELECTRIC CONSTANT 
SIGMA “ CONDUCTIVITY IN MHOS/METER 

REAP (lr*) EPSLN» STGMA 

*************************************** 


C 

c 

0017 
0018 
0019 
0020 
0021 
0027 
0023 
0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 50 

0037 

0033 

C 

C 

c 

0034 100 


COMPUTATION SECTION 

*************************************** 

SIGMA * SIGMA * l.E-U 

CHI - .053693 * D / LAMBDA**C1 

HI * HT<1) * 1.E-3 

H7 = HT<2) * 1.E-3 

DIS « (HI + RAD)**2 + <H2 Y RAD)**2 

DIS * SORT(DIS - 2#(HI + RAD) * (M2 Y RAD) * COS(D/RAD)) 

K = .002924 * LAMRDA**Cl 

K = K / SQRKSORK (EPSLN - l.)**7 + 3.6E75 * SIGMA**2 
1 * LAMBDA**?)) 

TE (NPOL.EG.t) GO TO 50 

K - K » SQRKEPSLN**? + 3.6E25 * SIGMA**? * LAMBDA**?) 

PSI * ATAN<EFSl N/(A.E12 * SIGMA * LAMBDA >> 

ESI * F'SI - .S * ATAN<(EPSLN - I.) / (6.El? * SIGMA * LAMBDA)) 
DELIA “ K * FXP(.I * ((3. * FI / .,) - PSI)) 

GO 10 100 

FSI « .5 * AIAN<(EPSLN - 1.) / (A.E12 * SIGMA * LAMBDA)) 

DELTA * K * EXPL.I * <<PI / .,) + PSD) 

FSI * PI / 2. - PSI 

RESIDUE SERIES 

K? * K**2 
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Table 1. (Cont.) 


0035 
OO 3 A 
003/ 
00 Ml 
00 TV 
0040 
004 I 
0041* 
0043 
0044 
0045 
004A 
004 7 
0040 
0049 
0050 
0051 

0052 

005 I 


0054 


0055 
005/* 
005 7 


0059 

00*0 200 
00*1 
00 A 2 
00 A 3 
006 4 
00A5 
00 A A 
00 A 7 
oo An 
0069 
00 70 
00 71 
0072 


00/3 

00/4 


00/A 
00/7 

00 7H 


OG/V 730 


K3 K**3 

K4 ■* K**4 
K5 - K#*5 

If <K.GT .O.A> 0(1 TO 200 

ANG1 - (PI / 4.) 4- PS I 

ANG? - 75. * PI / J BO, 4 3.* PSI 

ANG3 - 4. * PSI 

ANG4 = 75. * PT / 180. - 5. * PSI 
CANG1 -* COS(ANG1) 

SANGI - SIN(ANGI) 
gang? * cns<ANG2> 

SANG? - SINTANG?) 

MANS3 * COR < AN03) 

GANG3 =- SIN(ANG3> 

CANO4 - CfJR< ANG4 ) 

SAN04 - SlN <ANG4) 

TAUR(O) « . V7B 4 K*CANGl 4 1.237*K3*CANG2 .5*K4*CANG3 - 

l 2.755#K5*CANG4 

TALIR( 1) * 1.A22 4 K*CANG1 4 2. 1 63*M*CANG2 - ,5*K4*CANG3 - 
1 R,422*K5*CANG4 

TAUR(2) = 2.191 4 KfCANGl 4 2.921*K3*CANG2 - .5*K4*CANG3 - 
1 15.3A*K5*CAN04 

TAUR(3) “ 2.A937 4 KfCANGl 4 3.5919*K3*rANG? - .fi*K4*CANG3 
1 ?3.22AR*K5tCANG4 

'AUT(O) - 1.A07 - K*SANG1 - 1.237*K3*SANG? 4 .5*K4*SANG3 - 
1 2.755*K5*SANG4 

TAHITI) “ 2.810 - K*SANG.1 - 2.1A3*K3*SANG2 4 • 5*K4*SANG3 - 
1 B.422*K5*SANG4 

TAUT < ?) » 3.795 - KfSANG.1 - ?.921*K3*SANG? + .5*K4*SANG3 
1 15.3A*K5*SANG4 

TAU1(3) * 4.A A 3 3 - K*SANG1 - 3.5?t9*K3*SANG2 4 .5*K4*SANG3 
1 23.22AB*K5*RANG4 
GO rc» 300 

ANSI - 15. * FT / 1 BO * - PSI 
ANG? * 2. * PSI 

ANG3 “ 15. * PI /lBO 4 3. * PSI 
ANG4 - PI / 3. 4. * PSI 

FiCANGl = COS (ANSI.) 

RSANG1 •» SIN(ANGI) 

GANG? - COS<ANG?) 

SANG? “ SIN < ANG2) 

CAN03 * COS(ANG3) 

SANG3 « S J N< ANG3) 

CANG4 « COS < ANG4 ) 

SAN04 - SIN< ANG4) 

lAIJK(O) - .4043 4 . AI 8*8CANIi 1 /K * . 23A4*GANG;Vh? 

1 .0533*CANG3/K3 4 .002264CANG4/K4 

IAUR(l) - 1.280 t .194*HCANGl/k .00/J*SANGJ/K2 4 

I •0120♦C A NG 3/K3 - .00160*CANG4/K4 

I Alir < 0 > - .7003 . A183*HGANG1 /K 4 . 23A4*CANG2/K2 - 

1 .0533*SANG3/K3 - .00226*8ANG4/K4 

I AIM ( 1 ) = 2.232 ,194<>*BSANG1/K 1 .0073*CANG2/K? 4 

1 .0 I20*SANG3/K3 4 .00160*SANG4/K4 
Ml 230 I a 2.3 

1 AllR( 1 ) =-• I. 116*(I V .25>**C2 4 .2241*HCANG1 / 

1 < K* < I 4 . 25>**C2> 

TAU1(1) ~ 1•932#<I 4 .?5>**C? - .2241*KSANGl / 
l < k* < i 4 .>5 >**r:?> 

CUN IINUt 


00 HO 
0081 
0082 
0083 
0084 
0085 
OOBA 
008 7 
0088 
0089 


C 

C 

300 


399 


HUGH! GAIN FACT TIRO 


CflM ~ 50. * t AMttUA**C2 
SUN - 0. 

E * HAH. 3 * SilkT<CHI ) * 1.K3/UI0 
A < 0) “ .3582 * fc.XP<J * C2 * PI) 

A(l) * .3129 * EX P < - -. I * C! *PI> 

A(2) - .2903 * bXP<J * C2 * PI) 

A< 3) = .2/60 * FXP( - .1 * Cl *Pt> 

HO 399 N - 1»2 

CSO(N) =■ .036/4 * HI <N> / I.AMUHAMC2 
CONIINUb 
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Table 1. (Cont.) 


0090 no 899 I = 0#500 

0091 IF (l.Ul.J) GO TO 410 

0092 AX “ A(1) 

0094 TAG - CMPl X( TAUR< I) fTAIII (I ) ) 

0094 GO TO 500 

0095 410 If (K.GT.0.6) GO TO 450 

0096 TAURX = 1 ♦116*< If ,75>**C2 + K*CANG1 + I.4881*<It*75)**C2 

1 *N3*CANG2 - .S*K4*CANG3 - 3.9867*(1+.75)**C3*K5*CANG4 
009/ TAUIX = 1.932*<I+.75)**C2 - K*SANG1 - 1.4881*(I+.75)**C2 

1 *K3*SANG2 t .5*K4*SANG3 - 3.9867*<IF.75>**C3*K5*SANG4 
0098 GO TO 460 

0099 450 TAURX = 1.116*(1+.25)**C2 + .2241*DCANG1 / 

1 (K*<It .25>**C2) 

0100 TAUIX = 1,932*U + .25>**C2 -* .2241*DSANG1 / 

1 < K* <It* 25> **C2) 

0101 460 AX * .3440 * <-!)**<I+l) * EXP<-J * PI/3) / <I+.75)**C4 

0102 TAU * CMPLX< TAURX# TAUIX) 

0103 500 DO 700 N=1 #2 

0104 IF <HT<N).LT.COM) GO TO 550 

0105 IiIF CSCMN) - 2. * TAU 

0106 FI = - .1 * PI/4. + J * ( (SORT(IiIF) )**3) / 3. 

010/ F 2 = 1. - .1 * 5./<24. * (SORT < OIF) )**3) - 385./(1152. * 0IF**3) 

0108 F3 = 1. t J * 5./(24• * <SORT<OIF))**3) 

0109 IF <AOS<REAI <FD ) .GT.B9. ) GO TO 510 

0110 HUE<N) = F'2 * EXF‘<F 1 ) - F3 * EXP<-F1> 

0111 SSOJF * SURT<SQRT< OIF)) 

0112 HGF < N) = MGF < N) * AX / (DELIA * SS01F) 

0113 GO ID 530 

0114 510 HGF < N) = 1. 

0115 530 CONTINUE 

0116 GO TO 700 

0117 550 X - 4.E07 / LAMPDA 

0118 HGF < N) = t. t 6.283 * <1 / (DELTA * X**C1 ) - l./X) * 

1 HI <N> / I AM nil A 

0119 HGF<N> = HGF<N) - 39.48 * HT<N>**2 * <1. X**C2 * 

1 DELIA * TAU) / (DELTA * LAMBDA**2 * X**C3> 

0120 700 CONTINUE 

0171 SUMX * HGF <1) * HGF < 2) * FXP<J * TAU * CHI) 

1 / <7. * TAU - 1. / BELT A**2) 

0122 SUMX - SUMX * fc 

0123 SUM ~ SUM t SUMX 

0124 IF <1. fc>4 * APS < SUM >.GT.ABS < SUMX)) GO 10 900 

0125 899 CONTINUE 

0126 900 F ~ ADS < SUM) 

C ************************************** 

C 

C COMPUTE RECEIVED POWER IN DDU 

C 

r ************************************** 

012/ IF <1. AMBDA.l E .0. ) LAMBDA = 1. 

0178 IF (E.L T»0.) E « 1. 

0129 IF (PUR.IT.O.) PUR = 1. 

0130 RFCPUR = -158.513 t 20. * l OfilO<LAMBDA) + 20. * l.0G10<E) F 

1 10. * lOG10 < PUR) 

r: ************************************** 

c 

C OUTPUT SECTION 

C 

C ************************************** 

0131 SIGMA = SIGMA * 1.E11 

0132 UR1TE < 2 r1000) E r RECPWR # PUR r LAMBDA # D # HT #NPOL # EPSLN # SIGMA 

0133 1000 FORMAf </#lXr'OUTPUTS'.///#IX#'FIELD STRENGTH AT R * '#4X# 

1 IF'EV. 3# IX# 'MICROVOI TS/METER' r//#1X» 'POWER AT R 5 'r7X»OPF10.2 
7»6X f ' D B W' t///t IXt'*************************'*///»lXf'INPUTS' * 
3///»lX#'T PUUER * '#11 XfF9.3»6X#'KILOWATTS'# 

4//# 1X #'WAVELENG fH = '#8X»F9.3#6Xr'METERS'»//rIXr'DISTANCE 
5 TJR = ' »7X/FF).3r6X» 'KILOMETERS' #//# IX# 'HEIGHT OF T - '#7X 
6#f7.1#RX#'MF1 ERS'f//»1Xf'HEIGHT OF R ='r8X#F7.1#8X# 'METERS' t 
7//UX» 'POLARIZATION « ' » 1 OX f11 # T3B# ' 0 = VERTICAL# 1 = 

H HORIZONTAl 'f//fIX# EPSILON = '#14X#F5.2#//#1X#'SIGMA = '»17 
VX#1PE8*2»3Xf'MHOS/MEIER') 

0134 END 
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Figure 1. Comparison of ground wave r.m.s. field strengths calculated by 
the MRC ground wave program (x) and by H. Bremmer. 1 The 
ground constants apply to sea water, the polarization is verti¬ 
cal, the transmitter power is 1 kw, and both transmitter and 
receiver are on the ground. 
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Figure 3. Comparison of basic transmission loss in dB calculated by the 
MRC ground wave program (x) and by K. Norton. 12 The compu¬ 
tations are for average land (o = .005 mhos/meter, e = 15), 
vertical polarization and both transmitting and receiving 
antennas nine meters above the ground. 





































accomplished with an earlier version of the program using Bremmer's origi¬ 
nal conventions for transmitter power and distance between transmitter and 
receiver, not the CCIR conventions now used.) Figures 1 and 2 show that 
the MRC ground wave program matches Bremmer's computations almost per¬ 
fectly. The very small differences are probably due to the greater accur¬ 
acy possible with a computer. (All of Bremmer's results were computed by 

hand.) In contrast, the comparisons in Figure 3 between MRC and 
Norton^ reveal small although significant discrepancies. MRC's results 
are consistently lower than Norton's. These discrepancies probably arise 
because Norton also considers atmospheric refraction while we do not. 
Refraction is discussed further in Section 4. 

The MRC results also match the world standard CCIR results from 
the Geneva conference of 1974,but are somewhat lower than the CCIR 
results from the Kyoto conference of 1978. 14 This is not surprising 

because the 1974 results were based on Bremmer's equations while the 1978 

results were based on Bremmer's equations but also included refraction. A 
set of CCIR ground wave propagation curves is provided for reference in 
the appendix. 

The approach used by the MRC program is quite similar to that 
used by the NUCOM/BREM program.^ 

Using computations from the MRC ground wave program, the curves 
in Figures 4, 5, and 6 were generated. These curves apply to a 1 kW 
transmitter and to wide variety of circumstances from sea water to average 
land, from a frequency of 15 KHz to 10 GHz, from a range of 2 km to 2000 
km, and from a receiver height of 0 m to 1000 m. They demonstrate the 
versatility of the program while not attempting to cover all possible 
parametric variations. Figure 6 also shows that the program fails to 
compute correct values for a high receiver height and short distance 
between transmitter and receiver. The departures from a smooth curve at 
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Figure 4. Sample MRC results for ground wave propagation over sea water 
at various frequencies in MHz. Both the transmitter and 
receiver are assuned to be at zero elevation and both antennas 
are vertically polarized. 


















-40 



Figure 6. Similar to Figures 4 and 5 except the frequency here is fixed 
at 30 MHz and the height of the receiver (in m) is varied. The 
transmitter height is fixed at 10 m and the curves apply to 
average ground. 






receiver heights of 400 and 1000 m are errors. The program generally 
fails at very short separations between transmitter and receiver. In 
these cases, a 1ine-of-sight or geometric optical model is more appro¬ 
priate. 


Figures 7 and 8 show one application of the MRC ground wave 
program to a hypothetical ground wave system. The ground constants used 
are appropriate to average ground in Germany while the transmitter charac¬ 
teristics used are typical of many commercially available products. The 
maximim range is determined by the noise level at the receiver and not the 
receiver itself. Two curves are shown in each figure corresponding to a 
best and worst case noise. The actual situation will probably be some¬ 
where in between the two extremes. Figure 8 shows the considerable 
improvement in range that is possible with nominal improvements in trans¬ 
mitter characteristics. 
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Frequency 



Max Range (km) 10 dB S/N 

Figure 7. Maximum ground wave range for a hypothetical system in Germany. 
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Frequency 



Figure 8. Sane as Figure 7 except for Improved transmitter characteris¬ 
tics. 
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SECTION 4 
REFRACTION 


This section discusses the refraction of ground waves by a 
non-uniform atmosphere. Figure 9 shows that refraction is a significant 
effect which, under most but not all circumstances, improves ground wave 
propagation. 

The variation of refractive index as a function of altitude can 
be expressed in several ways. The following equations show two 
possibilities 


u(h) = 1 + A exp (-Bh) 

and 

u(r) = J I-n + y 


(19) 


( 20 ) 


where A and B as well as n and y are arbitrary constants chosen to suit 
the particular location involved, h is the height above mean sea level, 
r is the distance from the center of the earth, and a is the radius of the 
earth. Equation 19 comes from the CCIR* 4 while Equation 20 comes from 
Bremmer.l Although 19 is now considered the preferable form, 20 is 
adequate in the lower atmosphere. We choose to use Equation 20 because 
Bremmer's equations for refraction readily follow from this equation. 


Introducing the parameter a. 








too BO 300 _ „ 3S0!m 


Figure 9. Brenner's computation of the effect of atmospheric refraction 
on ground wave r.m.s. field strength.* It is interesting to 
note how the conposition of the air changes the refraction. 
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it can be shown* that the electric field E at a distance D and height h 2 
from a transmitter at height h t over soil with a parameter 6 (see Equation 
9) can be reduced to the electric field, E 0 , that would exist without 
refraction in the following way: 

E(D,h 1 ,h 2 ,6) = a 2/3 E 0 (Da 2 / 3 , h^ 1 / 3 , h 2 6 1 / 3 , 6a 1 / 3 ) (22) 

In other words. Equation 6 can be used to compute the electric field with 
refraction providing the modified parameters shown in Equation 22 are 
used. It is also necessary to use the actual wavelength existing at the 
earth's surface x 0 /u(a) and not the wavelength in vacuum. 

To be able to use Equation 22, it is necessary to evaluate a. 
This can be done using the following formula for the index of 
refraction:*** 4 


jj « 1 + 77.6 E. + .373 _£ 
T T 2 


(23) 


where 

T = absolute temperature (*K) 
p = atmospheric pressure in millibars 
e = water vapor pressure in millibars 

Differentiating 20 and 23 with respect to height yields. 


u'(r) = 


au(r) 


(24) 


and 


ii'(h) = -77.6*10- 6 -E. £L + 7 . 7 .-. 6 . xl0 ' . b . jg. - .746 — — + :ll 73 . — (25) 

T 2 dh T dh T 3 dh T 2 dh 


Equation 25 can be simplified by use of the relation 
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(26) 


dp de 
dh dh 


-10 gp 


where the unit of length is chosen to be 100 m and 

g = acceleration of gravity in cm/sec 2 
p = density of air in g/cm 2 . 

Substituting Equation 26 in 25 yields 


y ‘ (h) = - jZI p + (-77.6 xlO - 6 + ——) If!* (27) 

T T T dh 

- (77.6xl0“ 6 p + .746 i) — — 

T T 2 dh 

Because refraction takes place close to the earth's surface, p(r) and 
p'(r) can be adequately approximated by conditions at the earth's 
surface. Because 

Y-n = u 2 (a)-l ” .00058 « 0 (28) 

Y-n can be neglected. Then (24) can be written 

u‘(a) = - — (29) 

a 

Equations 21, 27, and 29 can be combined to provide the following solution 
for a in terms of atmospheric conditions near the earth's surface 


a 


1 - 



(77.6x10" 6 + -1^) 


1 de 
T dh 


( 30 ) 


+ (77.6xl0- 6 p + .746 -) — — ] 
T T 2 dh 
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Equation 30 when used in conjunction with Equation 22 is an adequate 
formal solution of the refraction problem which can be readily adapted for 
use on a computer. Using a better approximation than Equation 29 will 
lead to a slightly better solution of the refraction problem. 

Without a computer, the formal solution of Equations 22 and 30 
can still be used by employing several additional approximations. Using 
the approximation in Equation 28, Equation 22 can be rewritten 

ElD.hjjhj) = (1-n) 4 ^ 3 E o [n(l-n) 2 / 3 , h^l-n) 1 / 3 , h 2 (l-n) 1 ^ 3 ] (31) 

This equation further assumes that 6 (given in Equation 9) is small. That 
is, 

|6| = K « 1 . (32) 

This occurs for short wavelengths from VHF and shorter (perhaps also for 
portions of the HF band depending on the ground constants used). In any 
case. Equation 31 is considerably simpler than 22. 

Without refraction the electric field depends on D principly 
thru the parameter, x » given in Equation 7: 

X = (2ir) l/3 D/(a z/3 X l/ 3 ) = .053693 D^/^ 1 / 3 (33) 

From this expression, it is evident that changing D to D(l-n) 2 ^ 3 is 
equivalent to changing the radius of the earth to a/(l-n) and holding D 

constant. The effective radius of the earth can then be written 

a = a / (.766 - (.0681 + .00237 e) — + .306 —} (34) 

e L dh dh J 
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using Equations 27 and 29 and evaluating T and p at standard temperature 
and pressure (STP). 

Equation 34 is a most interesting expression. It is evident that 
for a dry atmosphere with an adiabatic temperature gradient (dT/dh «-l), 
a e « 1.20a. For saturated air, a e « 1.55a. For mean meteorological 
conditions a e *> (4/3)a, which is an often used value. Under some condi¬ 
tions a e < a, and the refraction will actually hinder ground wave propa¬ 
gation. In contrast, a temperature inversion (dt/dh > 0) will considerably 
enhance radio propagation. 

To compute the resulting electric field in the case of refrac¬ 
tion using the approximations above, it is only necessary to take a graph 
of E versus D without refraction and multiply E by (1-n)* 4 / 3 and D by 
l/(l-n) 2 / 3 . If h L and h 2 are non zero, then they must be similarly 
reduced as shown in Equation 31. The graph in Figure 9 was prepared by 
Reference 1 using this technique. 








SECTION 5 

PROPAGATION OVER IRREGULAR TERRAIN 


Up to this point, we have assumed that the earth is a smooth 
homogeneous sphere. Although the earth varies substantially from this 
idealization, the results expressed in Equation 6 are applicable to many 
situations. Nevertheless, it would be useful to compute field strengths 
for those cases where Equation 6 is clearly not valid. For instance, 
propagation over mountains or propagation from land to ocean is not well 
described by Equation 6. The literature on these subjects is extensive 
and well beyond the scope of this paper. We arbitrarily choose to discuss 
propagation over obstacles and to refer the reader to other work for 
discussions of propagation over inhomogeneous or inductive terrain. In 
all of these cases, the results occasionally differ from simple intuitive 
extensions of Equation 6 because of the unusual nature of diffraction 
processes. For instance, a large obstacle between a transmitter and 
receiver can actually increase the received signal strength. Such an 
effect is called "obstacle gain" and is explained below. 

Propagation over inhomogeneous terrain can be handled in various 
ways. The Suda method^ involves averaging the ground constants for var¬ 
ious segments of the terrain while the Millington method^ 7 *! 8 involves a 
geometrical mean of the electric fields due to the various segments. The 
Suda method is most appropriate where the ground*constants do not change 
substantially while the Millington method is applicable to substantial 
changes as long as the field is measured well away from the interface 
region. Many other references have considered these approaches as well as 
others. References 15 and 19 through 31 discuss the subject in detail. 
In addition. Reference 32 discusses propagation over inductive terrain. 
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Propagation over obstacles has also been covered extensively in 
the literature. The problem can be expressed in a formal manner in terms 
of a two dimensional integral equation of the following type:^»33 


W(p) 


1 


_1 

2tt 


// W(Q) e ik o( TQ + QP - PT) -EL 

QPxQT 


x [y(Q) - 



3QP 

3n 


(35) 


W is the ratio of the actual amplitude of the Hertzian vector to its value 
for the case of a perfectly conducting plane earth. The distances TQ, QP, 
etc. are indicated in Figure 10. The function y is related to the tilt of 
the' field as well as to so-called surface admittances. The integral 
extends over the entire irregular portion of the earth's surface except 
for an infinitesimal area around P. This equation can account for both 
irregular and inhomogeneous terrain. 


Equation 35 can be reduced to a one-dimensional equation 
provided that the irregular surface of the earth does not deviate too 
greatly from an average level plane or that the various ground parameters 
do not vary too greatly perpendicular to the line connecting the 
transmitter, T, and receiver, P. A saddlepoint approximation simplifies 
35 to 


Wj(x) 



J* w l ( . ^ [y i( S) + ik Q 

0 /c(x-e) 


sin 


<hU)] 


e ik 0 (TQ + QP-PT) dc 


(36) 


This expression shows that: 

a) Local changes of the ground constants and of the terrain 
profile have a similar effect on Wj. The term in square 
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brackets shows that y» representing ground constants, 
and a, representing terrain profile, have a similar 
effect. 


b) The terrain near the transmitter and receiver has a 
relatively greater influence due to the existence of the 
weighting factor l//£(x-c). 

It is possible to reduce Equation 35 to Equation 36 because of 
the dominating role of integration points, Q, near stationary values of an 
exponential phase factor. The region surrounding the area of a stationary 
phase is the first Fresnel zone. Figure 11 shows a drawing of this zone. 
Terrain irregularities and inhomogeneities produce significantly different 
results if they occur inside or outside this first Fresnel zone. Inside 
the zone, disturbances can prevent production of the main field at the 
receiver. This is the field that would have occurred had no disturbances 
been present. In contrast, disturbances outside the zone will not signif¬ 
icantly alter the main field but may add some perturbations. If the 
transmitter and/or receiver are elevated, the same reasoning applies but 
the first Fresnel zone is somewhat different. 

Major obstacles between the transmitter and receiver, such as 
the ridge shown in Figure 12, can produce a surprising result known as 
"obstacle gain". The ridge is assumed to be within the first Fresnel zone 
and well above the shortest path between the transmitter and receiver. In 
this case, the ridge acts approximately like the absorbing straight edge 
often considered in the theory of optical diffraction.2 The field at or 
beyond the ridge can be estimated from the Cornu spiral shown in Figure 
13. The starting point for the spiral is fixed at the lower point of the 
spiral, and the various chords, a through g, correspond to the magnK'de 
and phase of the electric field. Only the end point of the chord changes 
with x. (x is the coordinate perpendicular to the straight edge and 


41 









I 


parallel to the wave front.) In the geometrical shadow region (-<*> < x < 0) 
the length of the chord increases steadily as indicated by points a 
through d. d corresponds to the boundary of the shadow region. From this 
point, the chord keeps increasing in length until it reaches a maximim at 
f and then begins an oscillatory behaviour approaching the other end of 
the sprial. Figure 14 shows this effect explicitly. 

Behind the diffracting straight edge the field is given approxi¬ 
mately by 


- +- ( 6 /) 

cos(t2) cos(£H) 

where u represents the electric field when it is parallel to the edge and 
the minus sign applies. When the electric field is perpendicular to the 
edge, u represents the magnetic field and the plus sign applies. Equation 
37 assumes a cylindrical coordinate system with the z-axis parallel to the 
straight edge and the angular coordinate <{>=0 or 2n at the straight edge, 
a is the angle between the wave normal and straight edge. Because the 
value in square brackets decreases slowly with increasing <|>, the light or 
radio wave is diffracted far into the shadow region. Equation 37 is based 
on an approximation which fails at the shadow boundary <j> = a or <{> = -a. 

Equation 37 shows that the diffracted field beyond the ridge 
falls off as r -1 / 2 . This is less of a decrease than for ground wave 
propagation over a flat earth (r -2 ) or a curved earth (e _r//r °). In 
front of the ridge, the field decreases as r _1 which is also better than 
ground wave propagation. In other words, a sharp ridge can actually 

improve reception at a distant location behind the ridge . Because typical 
soil conditions produce considerable absorption of ground waves while a 
propagation path over a high ridge avoids the ground to a large extent, 
improved propagation is possible by utilizing a ridge. 


(1+i) g ikr 
4/nl<r' 
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Figure 13. 


Figure 14. 



Cornu spiral construction for determining the diffraction 
pattern produced by a straight edge. The axes labeled S and C 
represent the Fresnel integrals. 



Amplitude of the electric field at a straight edge. Zero 
corresponds to the boundary of the edge* negative values of w 
are behind the edge Mhile positive are in front of the edge. 
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Several computer programs have been developed to compute ground 
wave propagation over irregular and inhomogeneous terrain. These require 
a detailed knowledge of the terrain between the transmitter and receiver 
and involve long or short wavelength approximations. One program devel¬ 
oped by the Communications Research Centre 34 applies to frequencies at 
VHF or higher. Another, developed by ITS 33 applies to MF and lower fre¬ 
quencies but can be stretched into the HF. The basic difficulty with the 
HF is that the radio wavelength is close to the size of terrain features 
so that both long and short wavelength approximations have limited applic- 
abi1ity. 


Another difficulty in the HF band is the decided lack of experi¬ 
mental data to help verify computer routines. While ground wave data for 
rough terrain exists at MF and VHF, it is surprising that little exists at 
HF. Perhaps the HF has been ignored because it is used principly for 
longer distance sky wave propagation. In any case, data at HF is essen¬ 
tial to further our ability to accurately predict HF propagation over 
rough terrain. 

Several investigators have collected data over rough terrain at 
frequencies outside the HF. 33-4 ^ gf these the work by ITS is perhaps 
the most thorough. 35 * 33 This data has been used to show that the ITS 
program called WAGNER is indeed able to predict ground wave propagation 
over rough terrain below a few megahertz. 
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APPENDIX 


The following curves, considered to be the best available 
predictions of ground wave field strengths, are reproduced from the CCIR 
conference of 1978 in Kyoto, Japan. They come from the corrigendum to 
volume V, Propagation in Non-Ionized Media , and are dated February 25, 
1980. 


The curves refer to the following conditions: 


a) smooth homogeneous earth 

b) transmitter and receiver are on the earth 

c) transmitting antenna is an ideal Hertzian vertical 
electric dipole (nearly identical to a vertical antenna 
shorter than one quarter wavelength) 

d) transmitter power is defined in exactly the same way as 
indicated in Section 3 (basically 1 kw). 

e) refraction is accounted for by assuming an atmosphere in 
which the refractive index decreases exponentially with 
height. 

f) ionospheric reflections are excluded. 

g) distances are measured around the curved surface of the 
earth. 
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